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ABSTRACT 


The  work  under  this  contract  has  been  computational 
and  analytical  services  in  support  of  operations 
analysis,  applications  research  and  post  flight 
analysis  of  stratospheric  scientific  balloon  programs 
conducted  by  the  Air  Force  Cambridge  Research 
Laboratories  (AFCRL) .  Geophysical  and  flight  data, 
furnished  by  the  Government,  were  utilized  for 
transcription,  analysis,  graphing,  and  mathematical 
computations.  The  work  performed  includes  hand 
and  machine  plotting  and  the  writing  of  several 
computer  programs. 


COMPUTER  PROGRAMS  ANU  COMPUTER  RUNS 


1.1  A  printout  of  VHF  Omnidirectional  Range  (VOR) 
data  for  approximately  60  stations  using  a  program 
previously  written  at  Emmanuel  College  was  made. 

The  program  uses  geographic  and  magnetic  positions 
of  receiving  stations  to  calculate  latitude  and 
longitude  crossings  of  any  VOR  signal  received  from 
a  balloon.  This  work  was  accomplished  in  support 
of  an  in-house  research  effort  to  develop  balloon 
locating  systems  that  make  use  of  existing  naviga¬ 
tional  aids.  The  VOR  network  was  utilized  to 
provide  position  data. 

1.2  Programs  were  written  to  reproduce  and  analyze 
a  computer  tape  of  Weather  Bureau  Aerological 
Network  (WBAN)  data  for  kwajalein.  Wind  data  for 
selected  criteria  were  printed  out  using  the 
computer  program.  These  data  provided  meteorological 
input  to  feasibility  studies  for  a  Simulated  Large 
Explosive  Detonable  Gas  Experiment  (SLEDGE)  conducted 
for  the  Defense  Atomic  Support  Agency  (DASA)  by  AFCRL. 

1.3  Several  analyses  were  made  using  Rawinsonde 
tapes  supplied  by  AFCRL.  The  tapes  required  much 
editing  before  analysis  could  be  performed.  Some 
of  the  analyzed  data  were  plotted  by  an  electronic 
plotter.  This  work  provided  a  statistically  valid 
wind  data  base  to  be  used  in  balloon  trajectory 
analyses  for  balloon  programs  conducted  at  the 
White  Sands  Missile  Range.  Feasibility  studies  and 
operational  planning  data  were  provided  for  Project 
698AJ  recovery  tests,  "Study  of  Atmospheric  Aerosols" 
also  under  Project  698AJ,  and  the  National  Aero¬ 
nautics  and  Space  Administration  (NASA)  Cosmic  Ray 
Ionization  Spectrograph  Program  (CRISP). 

1.4  Full  printouts  were  made  of  data  for  balloon 
flights  from  the  AFCRL  launch  sites  at  Chico, 
California  and  Holloman  AFB,  New  Mexico  after 
compilation  and  analysis.  (Sec  Section  4).  This 
effort  used  a  program  previously  written  by  Emmanuel 
College.  Partial  printouts  of  selected  criteria 
were  made  for  balloon  reliability  studies. 


KEYPUNCHING 


Various  keypunching  done  included  FORTRAN  programs, 
meteorological  information  for  Standard  Atmospheric 
Tables,  and  balloon  flight  data. 


ANALYSES  AND  PLOTTING 


3.1  Wind  da'ta  were  extracted  from  microfilm  and 
analyzed  for  an  Advanced  Research  Projects  Agency 
(ARPA)  study. 

3.2  Winter  stratospheric  winds  at  Chico  and  Holloman 
(7  years  of  data)  were  analyzed  and  plotted  for  use 
in  flight  operational  planning. 

3.3  Surface  wind  durations  for  world-wide  geo¬ 
graphic  locations  were  analyzed. 

3.4  Time-altitude,  azimuth  angle  and  elevation 
angle  data  were  reduced  and  analyzed  and  trajectory 
plots  made  for  in-house  research  studies  on  balloon 
station  keeping  concepts. 

3.5  Wallops  Island,  Virginia,  launch  projections 
and  trajectory  dispersion  patterns  were  computed 

in  support  of  balloon  operations  conducted  by  AFCRL 
for  the  NASA  and  for  VOR  test  flight  planning. 

3.6  Stratospheric  wind  statistics  were  computed  for 
White  Sands,  Cape  Kennedy  and  Point  Mugu  as  input 
for  general  trajectory  studies  for  Air  Force  high 
altitude  balloon  programs. 

3.7  Vector  analysis  of  15-minute  wind  data  for  two 
balloon  flights  using  azimuth  angle  and  1IDO  were 
made  for  research  studies  of  minimum  wind  field 
fine  structure  in  the  lower  stratosphere. 

3.8  Tracking  analysis  of  two  balloon  flights 
consisted  of  checking  and  correcting,  if  necessary, 
VOR  readings  for  station  and  position;  determining 
approximate  location  of  the  balloon  using  triangu¬ 
lation;  plotting  position  using  VOR  information; 
and  comparing  results  with  Federal  Communications 
Commission  (FCC)  data.  Radar  data  were  also  plotted. 
These  data  provided  information  on  the  accuracy  of 
the  VOR  balloon  navigational  system. 

3.9  Miscellaneous  data  for  the  Balloon  Field  Test 
Director  Handbook  were  collected  and  prepared  for 
compilation  to  assist  the  Field  Test  Director  in 
carrying  out  balloon  operations. 

3.10  Balloon  Behavior  Study  (See  Section  6). 

3.11  Wind  Characteristics  Study  (Sec  Section  5). 
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3.12  M-hour  maximum  computations  were  made  from 
White  Sands  wind  data  for  use  in  tethered  balloon 
studies  and  operational  planning. 

3.13  balloon  gore  length  and  crease  distance  for 
various  types  of  balloons  were  computed  for  use  in 
balloon  design  research. 

3.14  Wind  data  at  10,000  feet  for  many  parts  of 
the  world  were  used  to  plot  m-hour  maximum  curves. 
These  curves  became  part  of  a  tethered  balloon 
study  for  the  Mid- Course  Surveillance  System  (MSS) 
p  rog  ram . 


4.  INVENTORY  LISTING  OF  FLIGHT  DATA  FOR  POST  FLIGHT 
ANALYSIS 


4.1  Information  was  extracted  from  post  flight 
summary  forms. 

4.2  Information  was  checked  for  consistency  and 
accuracy  with  data  in  flight  folders. 

4.3  Data  were  punched  on  IBM  cards  for  further 
processing  using  a  program  written  at  Emmanuel 
Col  lege. 

4.4  Multiple  copies  of  the  full  printout  were 
made  for  balloon  flights  conducted  at  Holloman 
(1962-68)  and  Chico  (1964-68). 

4.5  Selected  printouts  were  generated  to  analyze 
balloon  performance. 


5.  WIND  CHARACTERISTICS  STUDY 


5.1  Wind  velocities  were  extracted  from  rocketsonde 
data  (1959-67)  for  eight  stations  at  three  heights 
for  all  months.  Data  were  interpolated  and  units 
converted.  These  data  were  used  in  long  duration 
balloon  trajectory  studies  at  110,000  ft,  120,000  ft, 
and  130,000  ft. 

5.2  Statistical  analysis  of  wind  data  for  selected 
months  were  made  for  research  studies  of  trajectory 
dispersion  models. 
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5.3  Plotting  of  winds  for  several  altitudes  and  locations 
at  various  times  of  the  year  was  done  on  weather  maps 

to  determine  optimum  balloon  flight  scheduling. 

5.4  Coefficients  of  variation  for  speed  and  direction 
were  calculated  for  input  into  trajectory  studies. 


6.  BALLOON  BEHAVIOR  STUDY  -  RESEARCH  CONCERNING  TIIF.  VERTICAL 

floT  ToTT  "oF”b  aTl'oon'S -  - 


6.1  Baracoder  sheets  were  corrected  and  time-altitude 
data  generated  for  39  flights. 

6.2  Standard  Atmospheric  Pressures  for  various  altitudes 
were  calculated  from  tables. 

6.3  Temperatures  were  taken  from  WBAN  charts. 

6.4  Pressure  vs  temperature  charts  were  plotted  and 
corresponding  temperatures  from  pre-flight  radiosonde 
runs  were  determined. 

6.5  All  data  were  keypunched  for  analysis. 

6.6  Three-point  means  for  elapsed  time,  altitude,  and 
temperature  were  calculated  and  plotted  for  smoothing 
purposes . 


7.  BALLOON  FLIGHT  SIMULATION  PROGRAMS 


General:  There  are  four  programs  named  respectively 

SIMBALL,  BALLFLT ,  PARAFLIT,  and  MAGMED.  All  programs 
operate  on  input  wind  data  tape  (blocked  CC-808  described 
below  with  blocked  backup  tape)  to  be  mounted  on  tape  3. 

Input  parameters  for  all  four  programs  are  the  same 
except  as  noted  below  under  individual  program  description. 

TOPALT,  ISO,  IED--(F10. 0,214) 

TOPALT- -max imum  altitude  for  simulation  (feet); 

if  negative  indicates  flight  parameters 
are  changing  (see  below). 

ISD,  IED--start  date  and  end  date  for  simulation 
(inclusive);  first  two  digits  are  month, 
last  two  are  day. 
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If  TOPALT  is  negative,  the  following  are  expected: 

New  flight  parameters  (NAMHLIST  /ALTY/T I  ME , STOP) . 

New  TITLE  (13Ab)  for  printed  or  graphic  output. 

F 1 ight  parameters : 

TIME--Asccnt  rate  (feet  per  minute)  between 

(N-2)  thousand  feet  and  N-thousand  feet 
goes  into  TIME (N/2000- 1 ) .  Ascent  rate 
is  converted  into  time  in  minutes; 
alternatively,  ascent  rates  less  than 
50  fpm  are  assumed  to  be  time  in  minutes. 

STOP--Float  time  (minutes)  at  N-thousand  feet 
goes  into  STOP (N/2000- 1 ) . 


Programs  compute  statistics  and/or  solve  problems 
associated  with  balloon  flights.  Below  are  listed 
two  of  these  with  abbreviated  description  and 
idiosyncracies  if  any: 

7.1  SIMBALL--Basic  balloon  flight  simulator  described 

above.  Prepares  CRT  output  with  individual 
flights,  statistical  summaries,  50t  and  901 
probability  ellipses. 

Idiosyncracy :  Calls  PLOTID,  use  CRTPLT 
if  needed. 

7.2  MAGMEU--Program  to  compute  when  balluon  crosses 

point  86.4NM  west  of  launch  site. 

Idiosyncracies:  None. 


Following  are  listings  of  these  two  aforementioned 
prog  rams . 
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□ 


> 


1.  ID  *  5079  CANTO*  SJMRALL  *-l 

hSSTUP  3  CC-aOS* NORING 

6  I “J05 
1  I PFTC  MAIN 

VR< A«n,C)«(A-P*C/FN)/(FM-l #0) 

V  AR  ( A *  ? ) *VR!A»H*B) 

DIMENSION  CARO  180) 

DIMENSION  EZ 1 2  «  2 1 7 ) 

DIMENSION  TIMF  <  99) « STOP! 99) * AL  r ( 9  *>  »  VX !  99 ) « V V 1 99 ) 

NAMELIST  /ALTY/ TIME., TOP 
DIMENSION  TITLE!  13)  »  T  TL  1  2  > 

COMMON  T1TLC*  TTL 

COMMON  EX  «  E  V  «  SU”X  «  t »  *>  « I?  •  “  .  I  •  X  A  «  x*  .  I  i  «  I  EO  «  \f^RX  «  VARY  «Nr  L  I  Tb  *  h  N 
R  I  MEN:!  I  ON  TTT  (  P  *  .36 1  ) 

EOU  l  VALENCE  CTTT  <  1  « 362)  *2/.  ) 

DATA  TIME»iTOP/l»S*S#0«J#0«A#C«9av  .  *0 **9*3*0/ 

P  At)*37*2'.'b7Vc.»b 
I  0 A TE= 1 300 
'ALL  aLOT I  0 ( 1 .0  ) 

call  reread 

TEST  I  *  1 *0F7-1 .0 

1  RE  AO  (  *5  «  1  00  )  TOPALT  «  I bD ♦ 1 EO  « I  RE  W 

READ  ( 99  « 1 5 1 5 ) CARD 
WR l TF ! 6* 1 516) CARD 
IMS  FORMAT!  80A1) 

1516  FORMAT! 1X«*0A 1 ) 

WRITF!99*1 1A)  I  SO  t IED 

1  1  A  FORMAT! IA,AH  TO  IA) 

READ  199 1 1 15) TTL 
in-  FORMAT !F10»0«2IA«I1) 

IF!  ISD.LT.  IDATE  IRFVxlND  3 
IF ! TOPALT *GT *0*0 )G0  TO  68 
IF!T0PALT«EO#0.0)CALL  endplt 
READ ( 5  # ALTY ) 

READ <5. 1 151TITLE 
115  FORMAT! 13A6) 

DO  928  I  *  1 » 09 

IF! TIME!  I ) *GT *50*0) TIME!  I )  *2000.0/ T I Mb C I ) 

9?8  CONTINUE 

TOPALT*- TOPALT 
00  929  1*1*99 
J*2Q0 0*(  1  +  1  ) 

929  ViRITE  16*930)  J«TIME<  I  )  «  STOP!  I  ) 

93  FORMAT! I  1 0 *F 1 u* 8 »F 1 0*0 ) 

REMIND  3 
HR  NFLITEsO 
SUMX*0.0 
SUMY* 0*0 
SUMXX*0*0 
SUMYY=0*0 
SUMXY=0*0 

IF! I  DATE. EO* ISO > GO  TO  3 

2  REAC ( 3 ) S 1 TE  *  I YR  •  I DATF « I  HR »NALT » ( AL i  ( l  I  » VX 1  I >  « VY 1  I >  « 1 ■ 1 *NALT) 

1 F 1 S I TE *NE *T  EST 1 ) GO  TO  996 

IOATE= 1300 
GO  TO  995 

996  IF! IDATE.LT* I  SO) GO  TO  2 
IF!  I  DATE  *L(Z*  I  EO  )  GO  TO  3 
995  I F  1  NFL  I  TF.*GT  *  1  )  GO  TO  A 
WR1TEC6* 101 ) ISO* IED 

101  FORMAT! 2 1HONO  FLIGHT  PFTW:  EN  I6«5H  AND  16) 
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GO  TO  1 

3  IF ( ALT (NALT )  iLT «T  OP ALT ) GO  TO  2 

■•MITf  (6.102)  I  DATE  .  1  YR  .  (  ALT  (  1  >  .  VX  (  I  )  I  V  Y  (  1  )  .  I  a  1  «NALT) 

1  O'*  FORMAT!  I H 1 2  I  6/ ( 6F2C •  6  >  ) 

N'Fl  I  TE  =  -iFL  I  TF  +  1 
MM  A  X  a  TOP AL  T 
I VFC  =  1 
T=  3,0 
G=  0  •  0 
I  =6000 

666  J= 1/20-0-1 

I  VSC  =  1NTERPCNALT. ALT .  1-1000*  I VEC ) 

FF= I  -  1000 

T  =  T +FVFCT ( TIME, ALT, VX» IVtC.FF, J) 

5  =  S+FVEC'I  (  T  I  ME,  ALT.  VY  .  I  VcC.FF  .  J) 

FP=I 

I VEC  = INTERP ( NAL  T. ALT .  1  .  I  VtC ) 

■  =  G*F  Vf  CT(  STOP.  alt.  W  .  I  vEC.FF  .  J) 

TsT+FVECT ( STOP. ALT. VX. I VEC.FF . J) 

I =1+2000 

IF ( I .Lr.NKAXI GO  TO  666 
T=T/60.0 
S=S/60.3 
-Z<  1  . NFL  I TE ) =  T 
-  2  <  2 . NFL  I TE ) =S 
SUMX=BUMX+T 
SUMY  =  SUM Y+S 
SUMXX  =  SUMXX  +  T  *  T 
S'JMYY  =  SUMYY+E*S 
SUMXY  =  SUMXY+T *S 
GO  TO  2 
A  CONTINUE 
FN=NFLITE 

VAPX=VAR( SUMXX . SUMX ) 
v ARY= VAR ( S JMYY. SUMY ) 

R=VR ( SUMX Y. SUMX .SUMY ) 

PS  I  =  0 .5* AT AN2 <  R+R. VARX-VARY  >  *RAD 

XA=SORT( < VARX+VAPY) **2-A.0+ < VARX*VARY-R*R) ) 

XB  =  0  «5+ ( VARX  +  VARY-XA ) 

XA=O.E*( VARX+VARY+XA ) 

X A  =  SORT ( XA ) 

X -3  =  SORT  (  XR  ) 

R=R/SQRT ( VARX+VARY) 

SUMX  =  SUMX/FN 

SUMY  =  SUMY/FN 

EZ( 1 .NFLITE+I ) =0.0 

E  Z ( 2  «  NFL  I TE+ 1  )=0.C 

EZ  < 1  .NFL  I  TE+2 )=SUMX 

E Z ( 2 . NFL  I Tfc  +  2 ) =  SUMY 

CALL  FLOTR<EZ.TTT ) 

GO  TC  1 
END 

« ISFTC  DP° 

SUBROUTINE  PLOTR  (EZ.T) 

DIMENSION  EZ( 2.217) »T<?,361 ) 

DIMENSION  T I TLE ( 1 3 ) ♦ TTL ( 2 ) 

COMMON  TITLE .TTL 

COMMON  EX, EY, SUMX. SUMY. R. PS  I .XA.XB.  I SD , I ED . V ARX ♦ VARY . NFL  I T£ . Fn 
DIMENSION  GARB (6) 

DIMENSION  EL <2. 361) 

DATA  Bi_/1H  / 


-n  a  r 
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CALL  ELLlPS(EL.XA.Xb) 

CALL  ROT  (P£  I  *  SUMX.SUMY.EL  «T«2  •  12) 

WRITE (6. 1 15 >T 
WR I TE ( 6 » 1 15 )EZ 

W  R  I T  E  <  6 » 1 1 7 ) SUMX.  SUMY  »R .  PS  I .XA.XB 
17  FORMAT  <  //(  1P1  ?J-  10*2  )  ) 

1  S  FORMAT(//( 1?F 10.3)  ) 
r-:=5#  inflite+36?  ) 

I  I  I  =  1 

CALL  SCALE (T  . R «  0. N« 1 ♦ 20. 0 « XM I N »0X ) 

WRITE (6. 1 16)NFLITE»XM IMiOX 
16  FORMAT  < 1  HO. I  5  » 2F20 . 6 ) 

WRITE (6. 1 16)  1  II 
DO  1  1  =  1  « 13 

I  F ( T 1 TLE ( I ) .NE.BL >N  =  6*I 
CONT I NUE 

Ns-N 

CALL  FPAV.FI  O.o.?. 25) 

CALL  AX  I  SC C .0  « T ITLE.N  .6.0 .0.0. XT  I N . CX . 2C.0  > 

I  I  1=  I  I  1+1 
WRITEI6. 116)111 

C  ALL  AXIS(0«0*TTL«  12.U.C.W0.0.  XM 1 N  .  JX  .  20  •  0  ) 

I  I  la  I  I  1+1 
WRITE <6. 116)111 

CALL  LINK EZ «EZ(2* 1 ) «  NFL  I  TF. .  2  » -1 t 4 « XM I N «DX « XM I N . OX ) 

I  I  1= I  I  1+1 
WRITE  <  6. 1 16)  I  I  I 

CALL  L I NE ( EZ (  1 .NFLITE+I )  .EZC2.NFLI  i E+ 1  >  1 2 » 2 » U » 0 « XM I N » DX « XM IN « DX > 
111*111+1 
WRITEI6. 116)111 

CALL  L I NE ( T  . T ( 2 « 1 )  . 36 1 . 2 . l> . J ♦ XV  I N . OX . XM I N . OX ) 

111=111+1 
WRITE <6* 116)111 

CALL  LINE! T.T ( 2. 1  >  *  2 » 360 .0 « w « XM I N.DX » XM I N . OX ) 

I  I  1= I  I  1+1 
WRITE <6. 1 16) I  I  I 

CALL  LINcITC 1  .91) •TI2.91 )  . 2  *  360 «  0 ♦ 0  « XM I N  *  OX . XM IN*  DX ) 

I  I  1  =  I  I  I  +  i 
WRITEI6. 116)111 

CALL  ROT  (PS  I  «  SUMX  ♦  SLiM  Y .  EL  ♦  T  »  1  •  18  I 
111*111+1 
WAITS <6* 1 16) I  I  I 

CALL  LINEIT  . T < 2 . 1 ) « 36 1 . 2 . 0 » 0 « XM I N .DX t XM I N » OX ) 

I  I  I« I  I  1  +  1 
WRITEI6. 116)111 
V ARXa  SORT ( VARX ) 

V/  ARY  a  SORT  (  VARY  ) 

X=-XM IN/OX 
CALL  PLOT ( X  «  0  *  3 ) 

I  I  la  I  I  1  +  1 

WR I TE ( 6 . 116)111 
CALL  PLOT ( X . 8 • C .2 ) 

I  I  la  I  I  1  +  1 
WRITE (6* 1 16) I  I  I 
CALL  PLOT ( 0  *  X . 3 ) 

I  I  la  I  I  1  +  1 
WRITE  <  6. 1 16)  I  I  I 
CALL  PLOT( 8.0 .X .2 ) 

I  I  1= I  I  1  +  1 
WRITE (6. 116)111 
WRITE (99. 127) 


1  =  T  (  BOX  ) 

'.^ITC(  Q9  .1^1  )  >wr,X«  -oVY 

■■'‘fiAT  <  >c,  io:j  >  gar- 

I  I  “OR" A  T  (  ^iMXra  AR  =  C‘7.2»7X«6HY4AR  =  e-7.2) 

'  '  "OR"  AT  <  EAr ) 

ALL  jLk  {*,»— l»0«—3) 
i  I :  =  1 1 if  1 
.AIT*  (  *>»  1  16  )  I  I  1 

“all  ^y.’jcl<o*g«c*  l  ps.ga^* u .3*) 

•PITFI^iK  *  )  VA^Xt  v’ARV 

1  ?0R^AT(9«jr,vi-x  =  F7«  y.«JX*9HSIG'A-Y  =  P7»2) 

''All  ‘i.OTtv 
I  I  I  =  I i if i 
VR  I  TE  (  f>  •  1  1 6  )  I  II 
READ  <  ««.  102  )GAR- 

CALL  SV-'^OL  (  0  .0  «0.  1  25  .GAL'S-  ,0  .34) 

•‘ r i  te  <^9«i  34  )xa,xh 

104  ►'OR'-*  A  T  (  9M3  I  of>‘A  —  A  =  F*  7#  °  *  ?X  «  I  OMA--1  s  C>7*9J 

CALL  PLCT(C«-0,9-,-2.) 

1 1 1  =  1 1  If  1 

-JR  I  TEC  6*  116)  I  I  I 

RE  A  3 <  R9«  lO^lGAR- 

CALL  Y'HCL  (  0  «  3  .0*  1  P5*GAW'_  «0  «3A) 

CALL  ~LOT< C«-C*2b.-3) 

I  I  I  *  I  l  I  +  1 

’•'=9  I  Tt£  (6*116)111 

>(R  I  TE  (  °9  *  1  05  }  R.Pb  I  «  \FL  I  TE 

l.Ch  -OO'AAT  <  .3  Hf*  =  P7*At2x«SHPSI  =  FS  •  1  »  *:X  *  3HN  *  I4*6H 

READ  <  9Q  »  1C?  )  GAR.3 

CALL  i  vy-,CL(  0«r.  »C*  1 b5»GAkb.0.34) 

RETURN 
r  NO 

*!«FTC  ELL 

EU^ROuT  INF  ELL  I  PC  (  E  •  A  «B ) 

'I  V'V^ION  F(2«361  ) 

RAO* 3. 1 41 5926b JLH9/ 1«0*0 
-'O  !  1^1*76) 

X  =  FLCAT (  1-1  )*RAO 
F.  (  1  »  I  )  *A*COS(  X) 

1  “(7,1  )  sp*S I N ( X ) 

RETURN 

c-Nji-) 

Tl“eTC  RT 

^UaROuTIN-  ROT ( TMFT*XO«YO*EL*T»S) 

DIMENSION  EL(  ?.%  361  )  *T(2*36J  ) 

MANSION  XMAT(  ?«?)  *X(?») 

RAQ  =  3*  1  4  IE92 6E3U6Q/  1  fsO*0 
T  T  *  T  HE  T*R AD 
TT  =  -T  r 

X’MT(  1  .  1  )  =  COS(  TT) 

X^AT  <  1 .?>  =  3lN<TT) 

XYAT (  3.  1  )=-XVAT<  1  %?) 

X A T  (2.5  )=XMAT<  1.1) 

X( 1 ) *xo 
X( ? ) =YC 
10  1  J= 1  »  36 1 
70  1  I  =  1  *7 
'r  (  I  «  J)=X<  I  ) 

ro  ?  <-  l  »? 

^  t<  i .  j )  =  t <  i » j )  +  i#xv a  r  (  i  •  .<  >  *  -:l <  < « u ) 

1  continue 


1 1 
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RETURN 
END 
%  I  '"’FTC  FV 

^UNCTION  FVECT  (  T  «  A  «  X  «  I  «P  •  1 

O  I  WENS I ON  T(Q) .  A(R)  «X<P) 

Tfc’VP*  (  y.<  1  +  1  )-X  (  I  )  )/(  A  <  1  +  1  )  “  A  (  I  )  )  *  {  —A  (  1  )  >+X  (  I  ) 
FVECr  =  TL'-‘.P+T<<) 

RE  TUP'1 
“NO 

r  !  -FT C  I  NTP° 

-UNCTION  I  N  I';_RP  (  *  n  A  «  I  «  J  ) 

",r*“N?!Fig  A(QQ) 

I  F ( J.OT.N) J  =  N 
I  F ( J,LT.  i  ) J=1 
C  =  A  C  J  ) 

I  •  t  •  : •  1  )  uu  tj  l 

I  F  (  <  ,  _,T  «  1  )  GC  J 

L  (  J+l  ) 

I  F  (  <K  •  G£\  •  I  )  GO  TO  1 
J  =  J+1 
''-0  TO  2 

IF(J.  ‘£0*1)00  TO  1 
J=  J-l 
SO  TO  .2 
1  !MTERP=J 
RETURN 
-NO 

«^ATA 

10000# 0  0 7010715 

TALTY  T  I  W-* *99*800.0#  ST0P»9v*  0*0  s 
TO  10  <  FT  AT  80C  FPM.  SCALES  ARt  NM 
!Cn00#0  071 6070 1 

lOOQO.O  08010615 
10  '00.3  C 8 1  60b J 1 
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■iio  *  50  7  >  CANTOR  MAGMED  K-l 

iStToP  3  CC  —  S08  « NOR  I NG 

S  1  3  JOB 
SiaFTC  A  I  N 

VkIAc-iUsI  A---*C/FN  > / (FN-l *0) 

7AU<A»L})=\/R(A*:^«r>) 

jimensicn  caudibo) 

J I  PENSION  _2<  2*217) 

3  1  ML  NS  1  UN  T  I  :  (  99  )  «  SToP  (  ‘>9  )  »  ALT  (  9B  )  «  VX  (  99  )  t  V  Y  ( P  9  ) 

NAMELIST  /A_TY/TIV.e.  *STCP 
DIMENSION  TITLE  (  13)  *T  TL  (  2  ) 

CO -'.DON  TITL^«TTL 

C  JiViMQN  EX«EY«  NUVX«SU;  Y«n«PnI  «  XA  «  Xn  «  IsO*  I  feu*  V«RX*VARY#NFl  lTc.«r  N 
DIMENSION  TTT ( 2  • 36 1 ) 

EQU  I  VALENCE  (  T  TT  <  1  *362  )  «E2> 

DATA  T  I  ME  «  STOP/  i*0«2*0«3*0*4*0«  95*3  •  0  «  V  9*C  •  0/ 

RAD  =  57.  .-937793 

DATA  1 XX « 1 BL *XMAX/ t HX « 1H  #-5184*0/ 

I UATE= 1 300 
CALL  REREAD 
TEST  1  =  1 *0E7-1 .0 

1  READ ( o • lOOlTOPALT* ISO* IF.U«IREW 
10  FORMA T(F 1C*0«2I4* I  1) 

IF( ISD*LT* lOATfe )REWIND  3 
1 F ( TOP ALT • GT  *0*0) GO  TO  SB 
REAOtb.ALTY) 

READ ( 5  « 1 1 5 ) T I TlE 
115  FORMAT ( 1 3A6 ) 

WRITE <6. 1515) 

TMMrO.C 
DO  928  I  =1  *dcj 

f F ( J | M^| ^ | . GT  j  5  0 . 0 ) T  I  ME  ( I )«2000*0/T I  ME ( I  > 

CONTINUE 
TODALTs-TC^ALT 
DO  920  1=1*09 
J  *2000* (  1  +  1  ) 

029  WRITE (6* 930) J.TIMEI I ) «STOP< I ) 

93-*  F  DRV  AT  (  I  10«F10*3«"1Q*0) 

WRITE<6» l E 1 5 ) T I TLE 

ID  15  F0RMAT( IHl* 13A6/20H0DATE  T I V5  Y  ) 

REWIND  3 
c8  NFLITE=0 

F(  IDATfe.EQ.  ISD)G0  TO  3 

2  READ  (3  )SITE*  I  YR  «  I  DATE  «  lHR»NALT*  C  ALT  (  I  )  •  VX  (  I  )  « VY  (  I  )  «  I  ■  1  «NALT> 

IF (SI Tt*NE.TEST 1 )G0  TC  996 

I OATE= 1300 
30  TO  995 

A96  1F( IDATE*LT* ISC )GO  TO  2 

I F (  I  DATE  *LE • I  ED ) GO  TO  3 
r>95  IF  <NFL1  TE.GT*  1  )  GO  TO  4 
WRITE <6* 101 ) ISO* I  ED 

101  FORMAT ( 2 1 HONO  FLIGHT  rETwEEN  I6«5H  AND  16) 

GO  TO  1 

3  IF(ALT<NALT).LT*TOPALT)GO  TO  2 

C  WRITE (6* 102) I  DATE* IYR* (ALT(  I ) « VX (  I )  «VY! I  ) «  I*1«NALT> 

10->  FORMA  f(  I  HI  2  I  6/ (  6F2C.6)  ) 
nflite=nflitf+i 

NMAXr TOP ALT 

I VEC  =  1 
T  =  0  «  0 


•>  i 
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□ 


! 

i 

! 

I 


8 

m 

O 

M 

ft 

l 

\ 


s=o*o 

1  =6000 

666  JS1/20C0-1 

I  VEC  = INTERP(NALT»ALT« I- 1 000 ♦ 1 VEC) 
p-F=  1  -  1000 

T  =  T+F  ,/ECT  (  T  1  ME  *  ALT  t  VX  «  !  Vi.  C  »FF  %  J  ) 

:d=S+F7ECT<  T  I  ME  •  ALT  «  VY  «  !Vc.C|FFiJ) 

i out  = i^l 

2F< T.GT.XMiXJOO  TO  667 
I  O'.JT  =  I  XX 
GO  TO  660 
66’’  CONTINUE 

1 i 1 jOw 

1 F (  I «  LF • NMAX ) GO  TO  666 
FF=TQPAL7 
1  sNMAX 

I  VEC  =  INTEPP(NALT«AL  r «  1  «  IVEC) 

X*FUNC ( STOP. alt . VX« I VEC«F F. J) 

STOP  <  J)  =  ( XMAX-T )/X 
TMI =TMM+STOP{ J) 

Y  =  S  +  FVF.CT  (  i> TOP «  ALT  t  VY  •  1 VEC«FF  »u) 

Y=Y/6  > • 0 

1  F  {  X  •  G*  0  )  Y=TEbT  1 
660  CONTINUE 

I VEC»I VEC+l 

WP!TE(6» 1 14 ) IOATE* I YRtTMl « 10UT  *Y 
114  COQMATI  1 5  «  I  2  «  F7*  1  *  A  1  «F7#  1  •  EX  •  I  3  «2  1 3 1  1 P7E  1  3  •  3 ) 
STOP( J)«0.0 

GO  TO  2 
4  CONTINUE 
GO  TO  1 
ENO 

Z  I  OFTC.  Fv 

FUNCTION  FVECTt T«A«X« I fP«K> 

01  MENS  I  ON  T(9).AI9)»X(9> 

TEVO=<X<  l  +  !)-X(I))/(A(I  +  l)-A(m#(=-A(  DKXI  I) 
FVECT  =  TE  MP*T  ( '<  ) 

RETURN 

end 

vIBFTC  INTTT 

FUNCTION  FUNC<T»A«X« l •PfK) 

01  MENS  I  ON  T(9)«A(9)«X<9> 

TEVIP=  <  XI  14-1  )-X<  J  )  )/<  A(  1  +  1  >-A<  I  )  )  *  ( P  — A  (  m+X(  I  ) 
pUN'C  =  TEM° 

RETURN 

END 

3I0FTC  IN-PP 

FUNCTION  I M7EWP <  N«  A  « 1 « J ) 

O I  MENS  I  ON  4(9P> 

IF( J.GT.N) J=N 
IF{ J.LT. 1 ) J=1 
*  <  =  A  {  J  ) 

IFIK.FQ,  I  )C,0  TO  1 
IF ( <  «  GT • I ) GO  TO  J 
t<<  =  A  {  Jt- 1  > 

I F ( <K • G£ • I ) GO  TO  1 

JS  J+l 

GO  TO  ?- 

.7  IF  (  J  »  r"  Q  •  i  )  GO  TO  1 
J  =  J-  1 
GO  TO  2 


1 


I NTERP= J 

3£TU»N 

rNT 


IS 


□ 


t,  >a'rA 

11  0  7.1  13831 

‘  AL TV  Tr.'P:n#Oi??*aOOiO  ?. 

FDV  TO  11CCOO  ft.  PROFILE  A 

i  i  "i  'o,o  o  v"  i  zpj  i 

i  Ai_  TY  T  I  vE-0»J»i?3*8uG«G«  30*600*0$ 

F.T«  xo  50  <-ft,  600  FPM  TO  lie  <-FV«  PROF  ILF  H 
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! 

I 

8 

m 

sJ 

« 

ft 
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PROJECTS  SUPPORTED  UNDER  THIS  CONTRACT 


1.  Air  Force 

6665 

8292 

698AJ  (Recovery) 
698AJ  ( AOHP) 

2.  DASA  -  SLEDGE 

3.  NASA 

Goddard  -  BAPE 
MSC  -  CRISP 
Langley  -  Viking 


4.  ARPA  -  1366 
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